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ABSTRACT 

A model is presented in which the strongest radio-emitting jet outflows are produced in black hole 
systems when the accretion is a geometrically thick {H/R ~ 1) inflow {e.g., ADAF, CDAF) and if the 
black hole is rotating. For galactic black hole candidates, the model naturally accounts for the observed 
correlation of jet outflow with the black hole hard emission state and predicts an association of strong 
jets with rapid black hole rotation. When extended to the supermassive case, the model accounts for 
the highest radio galaxy and quasar jet powers and provides additional theoretical support for the "spin 
paradigm" , which asserts that radio loud quasars are produced by Kerr holes and radio quiet ones by 
Schwarzschild holes. In some cases, the angular momentum and energy outflow in the jet may be large 
enough to significantly alter the structure of the accretion flow from that predicted by current models. 

Subject headings: black hole physics — galaxies: jets — galaxies: nuclei — hydrodynamics — MHD — 
quasars: general — radio continuum: galaxies — relativity 



1. INTRODUCTION 

There is now substantial evidence that, in binary black 
hole candidate X-ray sources (BHCs), a radio-emitting, jet 
outflow is produced only when the X-rays are in the hard 
state. That is, the jet is "quenched" when the dominant 
accretion disk emission is from soft (2-10 keV) photons 
and little or no emission is seen in the hard (20-100 keV) 
region. Sources which display th is behavior are GX 339-4 
flFcnder et al. 1999[ Render 19991), GRO J1655-40 (a t least 



in the first few 1994 outbursts; 
1915+105 (IHarmon et al. 1997 



X-3 ([VlcCollough et al. 
et al. 19991 ). 



I Harmon et al. 1995|) GRS 
|Fender et al. 19991), Cyg 



1999| ), and Cyg X-l flBrocksopp 



In addition, there is some evidence that those black hole 
systems which produce a strong jet contain Kerr (rapidly- 
rotating) black holes, although this result requires addi- 
tional confirmation. Recent attempts to measure the an- 
gular momentum parameter j = J / (GM 2 / c) {i.e., "a/M" 
in geometric units) indicate large values in the superlumi- 
nal sources GRS 1915+105 and GRO J1655-40 (j = 0.95) 
and smaller values in the non-superlum inal systems Cyg 
X-l and GS 1124-68 (|Cui et al. 1998|). Furthermore, in 



the supermassive black hole case this "spin paradigm" has 
useful advantages for explaining the difference between ra- 
dio loud quasi-stellar radio sources (QSRs) and radio quiet 
quasi-stellar objects (QSOs) and why giant radio sources 



occur only in elliptical galaxies ( Wilson fe Colbert 1995 
Meier 1999|) 



This letter shows that these properties are a natural 
consequence of black hole accretion and can be explained 
by combining three current areas of astrophysical theory 
into a single model: 1) modern disk accretion theory, in- 
cluding geometrically thick inflows at low accretion rate, 
standard disks at intermediate rate, and thermally unsta- 
ble disks at high rates; 2) standard jet-production theory, 
in which the jet power is a strong function of the poloidal 
magnetic field strength and of the angular velocity of that 
field; and 3) dynamo theory for the accretion viscosity, in 



which the poloidal field strength is a strong function of the 
disk height. Thin disk flows, despite their rapid Keplerian 
rotation, do not have strong vertical fields and, therefore, 
should not produce strong jets. Geometrically thick accre- 
tion flows, on the other hand, are expected to have strong 
vertical fields, but suffer from slower rotation. Only when 
the disk is thick and the geometry also rotates rapidly, as 
it does near a Kerr hole, do we obtain both strong fields 
and rapid rotation and, therefore, strong jets. 

2. SUPERMASSIVE BLACK HOLES — RADIO GALAXIES 
AND QUASARS 

2.1. Maximum Total Observed and Predicted Jet Power 

The maximum observed total jet power possible from su- 
permassive black hole systems will be used to constrain the 
models for jet produc tion. To do so, the upper envelope of 
the OLWL diagram (Owen fc Laing 1989) ; pwen fc White 
1991; Lcdlow & Owen 1996) for powerful radio sources is 
used to give a maximum radio power as a function of the 
host galaxy's red optical magnitude. This relation is given 
approximately by 



logP™r « -0.56M* 

with the radio power P in WHz -1 . 
mendy fc Richstonc corre lation for red magnitude (Kor- 
mendy & Richstone 1995) 



14.9 

Secondly, the Kor- 



log m 9 



-0.5M R - 11.3 



is used to determine radio power as a function of black 
hole mass, where nig = M,/10 9 Mq. Finally, a factor of 
- lO" 11 - 10" 12 Hz _1 flBicknell 1995| ), corresponding to 
a fraction e ~ 10~ 2 — lO -13 of the jet power emitted in a 
10 9 Hz bandwidth, is used to convert radio power to total 
jet power, yielding the following expression for an estimate 
of maximum total jet power 



j max 
L jet 



10 



46 , 



i 

2.5 



(1) 



ergs ulq 

Any mechanism of jet production must be able to achieve 
this maximum power, not only instantaneously but also 
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sustained over long periods of time, as well as scale roughly 
linearly with black hole mass. 

To compare this observed power with theoretical es- 
timates, jets from both standard disks and advection- 
dominated accretion flows (ADAFs) will be considered, 
comparing the results for both Schwarzschild and Kerr 
black holes. This discu ssion is similar to that of ( Livio 



Ogilv e~fc Pringle 1999| ; LOP), except that the jet power 



from non-spinning and spinning black hole systems will be 
contrasted, rather than comparing the power from differ- 
ent components (disk and hole) within a single system; and 
geometrically thick accretion flows will be considered as 
well as thin disks. As LOP point out, when the hole is spin- 
ning, the power generated by the mag netic field threading 
the disk (the Blandford-Payne process; Blandford & Payne 
19820and that from the fie ld threading the black hole (the 
Blandford-Znajek process; Blandford & Znajek 1977) will 
be comparab le, at best. However, as pointed out earlier 



(Meier 1999), even if the BZ process is neglected entirely, 
the jet power contributed by the field threading the disk 
alone still will be a function of the black hole spin, since 
the rotating metric contributes to the rotation of the mag- 
netic field. That is, as viewed from an outside observer at 
infinity in the Boyer-Lindquist reference frame, the disk 
angular velocity is a sum of its angular velocity relative to 
the local metric f2' (i.e., in the Zero Angular Momentum 
Observer or ZAMO frame) plus the angular velocity of the 
metric itself in the Boyer-Lindquist frame lu = —g^t/gipip 

fi = fl' + LU 

The structure equations for accretion flows (e.q 



Shakura fc Sunyaev 1973 ; |Novikov fc Thornc 1973 



Narayan fe Yi 1995) present a solution for the disk mag- 



netic field which generally is identified with the domi- 
nant azimuthal component B^_ However, models and 



simulations of jet production (Blandford fc Znajek 1977 



Blandford fc Payne 1982; Meier et al. 1997) show that 



it is the strength of the poloidal field component (B p = 
[B 2 R + B 2 ?] 1 / 2 ) that is largely responsible for the produc- 
tion of a jet, with the output power given by 

L jet = B* B*n 2 /4c (2) 

where Rq is the characteristic size of the jet-formation re- 
gion, and B p o and are, respectively, the poloidal field 
and angular velocity of the field in that region as measured 
by an external observer. Now, LOP have arg ued, and 
dynamo simulations tentatively have confirmed ( Khannj 
19990 that the strength of the poloidal magnetic field is 
limited by the vertical extent of the turbulent eddies. That 
is, Bpo/B<p w (H /Ro) n , where H is the vertical half- 
thickness of the disk at radius Rq, and the exponent n ~ 1, 
yielding an estimate of the jet power of 



J jet 



B 2 H 2 R 2 Cl 2 /4c 



(3) 



which can be used to calculate the instantaneous jet power 
once B,f,o, Hq, and flo are determined using a particular 
disk model. 

2.2. Jet Power from Standard Thin Accretion Disks 

First, standard disks will be considered. For illustrative 
purposes, the "middle region" solution will be used, where 



gas pressure and electron-scattering opacity are dominant. 
For the Schwarzschild case (u> = 0), a characteristic ra- 
dius of Rq = 7GM/c 2 (just outside the last stable orbit) 
will be used, along with the Keplerian angular velocity of 
Q' = l^/ 2 (GM/c?y 1 . For the Kerr case, since the disk 
extends inside the ergosphere (inside R = 2GM/c 2 ), a 
characteristic radius of X.hGMjc 2 will be used, along with 
luq = lo(Rq) w 0.6(GM/c 3 )~ 1 . This yields 



^Schw-SS 
■'jet 



10 42 - 6 ergs" 1 (c^fr 1/10 ml /W m 6 ^ 5 

rKerr-SS ^ 1Q 43.5 ergs -l (a £f J-l/10 ^J™ ^6/6 



(l + l.lj + 0.29/) 



(4) 



(5) 



where a^f is the viscosity parameter for standard thin 
disks in units of 0.01, m_i = M /Q.lMEdd is the accretion 
rate scaled to that producing one-tenth of an Eddington 
luminosity, and the Eddington accretion rate is defined 
as MEdd = 22MQyr -1 mg (assuming 10% accretion effi- 
ciency). One-tenth the Eddington accretion rate is used 
as the scale value since it is roughly the rate at which 
the low/hard to high/soft black hole state transition is be - 
lieved to occur (Narayan, Mahadevan, & Quataert 1998). 

As LOP point out, the black hole spin makes little dif- 
ference for standard thin disks. Within the Kerr system, 
it contributes only about a third of the power. And the 
ratio of power from a Kerr to a Schwarzschild hole is ~ 20, 
with most of this due to the deeper potential well in the 
Kerr case. Furthermore, the jet power for a 10 9 Mq hole 
is limited to ~ 7 x 10 43 erg s _1 — two orders of magnitude 
below our target in equation (|l|). So, while Shakura & 
Sunyaev disks might account for some of the weakest ra- 
dio sources (in particular, the "radio quiet" QSOs, when 
there is significant optical emission from the disk), they 
cannot account for the most powerful radio galaxies and 
quasars. 

2.3. Jet Power from Advection- Dominated Accretion 

Flows 

Next, ADAFs will be considered, as an example of a ge- 
ometrically thick accretion flow. Note that this case will 
include coronae of thin disks in which a dominant fraction 
of the gravitational binding energy of the accreting ma- 
terial is dissipated in the corona rather than in the disk. 
For these flows H ~ R 0l so B v0 ~ B M . With the ADAF 



solutions (Narayan, Mahadevan, & Quataert 1998), both 



the Schwarzschild and Kerr cases yield 



L 



-Schw-AD 
J jet 

Ken — AD 



j<'t 



where / 



lO^ergs- 1 ^^)- 1 ™^!/ 2 
10" 



46.0 -1 



( a -?/2) 1 m g rh- 1 g 2 



(0.55/ 2 + 1.5/j+j 2 ) 



(6) 



(7) 



is the ratio of the actual angular veloc- 



ity to that calculated by Narayan and Yi, and g = g *° Y 
is a similar ratio for the azimuthal magnetic field. Note 
that, as is customary, the viscosity parameter for ADAFs 
is now scaled to ~ 0.3 instead of 0.01. While considerably 
larger than the jet powers calculated for thin disks, these 
expressions are not simply (Rq/Hq) 2 times equations (Q) 
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and (0) . This is because ADAF azimuthal fields are weaker 
than in thin disks, and the rotation rate of the inflowing 
plasma is somewhat smaller (0.4 Keplerian, as determined 
by Narayan & Yi). In fact, f2 is very sensitive to the 
exact value of the internal adiabatic index: f could be 



jets, and the association of the onset of jet production in 
these systems with a hard emission state, is very similar 
in character to the need for both spin and advective accre- 
tion flow to produce extragalactic jets discussed above. In 
fact, diag rams drawn depicting this association (Fender e\ 



consic crably less than unity (r = 5/3 yields f = 0!), so al. 1999 ) can be qualitatively, and to a great extent quan- 



the results have been scaled to f z . Equation (||), there 
fore, is probably an upper limit for the jet power from an 
advection-dominated flow into a Schwarzschild hole. 

On the other hand, in the Kerr case, the actual values of 
B^o may be higher than that calculated by Narayan & Yi, 
since these authors did not take the field- enhancing s hear 
in the Kerr metric into account (see e.g., Meier 1999| ). It 
is not unreasonable to expect an enhancement of at least 
g ~ oj(Rq)/Qqn-y ~ 2.3, with a corresponding factor of 
5 or more increase in jet power. This alone would put the 
Kerr- ADAF jet power for a 10% Eddington accretion rate 
well above the target of 10 46 ergs -1 . Equation (f7j), there- 
fore, is likely to be a lower limit on the jet power from 
an advection-dominated flow into a Kerr hole. Unlike the 
thin disk case, then, the black hole spin in the ADAF case 
is very important, especially if / — > 0. The ratio of Kerr 
to Schwarzschild jet power (~ 20<7 2 // 2 ) is potentially 100 
or more. 

The jet power in equation ([?]) is a substantial fraction 
of the accretion luminosity 

r Kerr — AD 



= 0.80 {at f/a)" 1 (0.14/ 2 



""i/./ • j 2 )g 2 

(8) 

Consequently, the angular momentum and energy carried 
away from the ADAF will have a substantial effect on the 
structure of the inflow, which must be taken into account 
in the ADAF solutions. (This is in addition to any ther- 
mal wind that may be part of an Advection-Dominated In- 
flow/Outflow Solution o r ADIOS which also must be t aken 



into account. See, e.g., Blandford fc Begelman 1999 .) In- 
deed, the ratio could easily exceed unity, removing more 
angular momentum and rotational energy from the inflow 
than it has in the ZAMO frame. This could reverse the 
spin of the ADAF in the ZAMO frame, giving it nega- 
tive angular momentum, and ultimately spinning down the 
hole slightly when the accreting matter finally enters the 
horizon. Such a scenario is quite plausible in Kerr geom- 
etry and, indeed, has been seen in magnetohydrodynamic 
(MHD) num erical simulations of jet production near ro- 
tating holes (Koide et al. 200C). This process is a type of 
Penrose mechanism, with the disk particles scattered into 
negative angular momentum orbits by a torsional MHD 
wave (the jet) that, in turn, is "scattered" upward and 
outward to infinity. Similar in character to the Punsly - 
Coroniti "ergospheric winds" ( Punsly fc Coroniti 1990| ), 
this process can extract angular momentum and energy 
from the black hole even without the magnetic field thread- 
ing the horizon. 

To summarize, the most powerful extragalactic jet 
sources require both a geometrically thick inflow and a 
rapidly spinning black hole to reach the observed maxi- 
mum jet powers. 

3. STELLAR MASS BLACK HOLES — MICROQUASARS 

The possible association of black holes having poten- 
tially high spin with binary X-ray sources that produce 



titatively, understood using the above equations. For the 
disk structure in different accretion r egimes, the discus- 
sion below will adopt t he solutions of Che_n et al. (1995)1 
and |Esin et al. (1998)[ Equations (@) - (@) will be used 
to calculate the jet powers, with the appropriate scaling 
of mass and accretion rate. Note that the jet power pre- 
dicted by these equations is primarily kinetic energy and 
Poynting flux; the actual radio luminosity emitted will be 
much less. For example, during the first 1994 outburst 
in GRO J 1655-40, which was thought to have reached ac- 
cretion rates appr oaching in ~ 1, the total jet power is 

-|38-39 



10 3S 3y ergs 1 (Meier 1996|), while the radio luminosity 



is only ~ 10 32 ergs 1 (Harmon et al. 1995). 

For low accretion rates, to < m A w 0.1 (a^^, 2 ) 2 ac- 
creting black holes are in the low/hard state — steady 
accretion with a substantial fraction of the X-rays in the 
20-100 keV band. This corresponds to an ADAF solution, 
which is favored over a thin disk solution in this regime. 
As the ADAF is stable, and expected to produce a jet 
power of up to ~ 5 x 10 38 ergs" 1 m\ rh-i j 2 (equation 0), 
a continuous jet outflow that is roughly proportional to the 
accretion rate is expected. This is, in fact, seen in several 
microquasars ( Fender 1999 ). Note that, although propor- 
tional to the accretion rate, the power in this model comes 
mostly from the rotation of the black hole and not from 
the accretion energy of the infalling matter. However, the 
efficiency of extracting the rotational energy of the hole is 
proportional to to. 

For to > to a the disk enters the high/soft state 
( Narayan, Mahadevan, fc Quataert 1998 ), in which the 
hot ADAF largely disappears and the disk takes on an 
optically t hick, cooler structure sim ilar in nature to the 
standard ( ghakura fc Sunyaev 1973 ) thin disk solutions. 
Here the cool disk is modeled with the Shakura & Sun- 
yaev "middle" region. The expected jet power from such 
a thin disk with a weak poloidal field (equation ||) is 

~ 5 x 10 36 ergs" 1 to^ 10 to 6 ^ 5 — a hundred times weaker 
than in the ADAF case. This comp ares well to the jet 
suppression of > 35 seen in GX 339-4 ( Fender et al. 1999 ) 
when it goes into the soft state, and is some of the best 
evidence that a strong jet cannot be produced unless the 
accretion flow is geometrically thick and hot. 

For high accretion rates (to ;> 0.5 — 1.0) the disk enters 
the "very high state" , in which the disk is very unstable. In 
GRS 1915+105, for example, the disk goes through cycles 
of ~ 20 minutes in duration, which begin with a dramatic 
dip in the soft X-rays, followed by an increase in the hard 
X-ray emission for several minutes. Just as the hard emis- 
sion abruptly disappears, a radio jet is ejected, followed by 
a softening of the spectrum again. A possible explanation 
for this behavior is the onset of a thermal "limit cycle" , in 
which the disk oscillates between two unstable accretion 
states (standa rd soft and hard ADAF) . Deta iled modeling 
of this case by Szuszkiewicz fc Miller (1998) suggests that 



the disk first bloats to a time-dependent, hot, gas- pressure- 
dominated ADAF-like state; then it drains rapidly into the 
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hole; and finally a thin radiation-pressure-dominated thin 
disk re-fills, only to become unstable and bloat again. For 
a 10 Mq hole, the period of the limit cycle is on the order 
of ~ 1000 seconds, with the active ADAF phase lasting 
only ~ 1% of that time. 

There is some uncertaint y in this model, however. Stan- 
dard disk stability theory ( Shakura fc Sunyaev 1976| ) pre- 
dicts that unstable behavior should occur above the rather 
low accretion rate of m J> 0.02(a^.f mi)" 1 / 8 , when the in- 
nermost part of the disk becomes r adiation pressure dom- 
inated. Indeed, the unstable disk of 3zuszkiewicz & Millci 
(1998| had m = 0.06. This rate is so low, however, that 
it predicts that all disks with m > rriA should be unsta- 
ble; that is, there should be no stable high/soft state even 
in the case of X-ray binaries; at most the disk should be 
in a cool state for a portion of the limit cycle (20 min- 
utes). This is in direct conflict with the observations: the 
disk stays in the high /soft state often for months. Using a 
more careful analysis, Chen et al. (1995) derive a modified 
transition accretion rate of 



rh R = 0.3(a£fmi) 



(9) 



(minimum of the solid lines in their figure 2) , above which 
limit cycle behavior should set in. Equation (|) agrees 
well with the observations of GRS 1915+105 and GX 339- 
4. Note also that the solutions of Chen et al. (1995) also 



predict that unstable behavior should continue for accre- 
tion rates well in excess of the Eddington rate (m <J 20 or 
so), above which a steady, radiation-pressure-dominated 
ADAF should set in. 

Ascribing a high accretion rate (m ~ 1) to the hard 
temporary ADAF state, equation (m predicts a short- 
term jet power during the very high state of ~ 6 x 
10 39 ergs -1 to 3 /2 mg 2 for a rapidly-rotating black hole. 
Again, allowing for possible metric-shear enhancement of 
the magnetic field, instantaneous jet powers well in ex- 
cess of 10 40 ergs _1 in GRS 1915+105 are possible in this 
model. 

4. A GRAND UNIFIED SCHEME FOR ACTIVE GALACTIC 

NUCLEI 

Finally, we turn once again to the supermassive case. 
The above concepts allow us to begin to build a scheme 
that unifies radio loud and radio qu iet AGN. It is similar to 
that discussed earlier ( Meier 1999[ ), but with the addition 
of geometrically thick advection-dominated flows. 

Since radiation pressure is even more dominant in super- 
massive black hole disks than in stellar mass ones (rhu oc 

m 1 1 ), it is likely that ttlr will be less than for some 
large hole mass (from equation 0, for M. > 7 x 10 4 Mq). 
That is, for supermassive holes, a stable high/soft state 
should not exist. Unlike before in the stellar mass case, 
the possible lack of a high/soft state is not in conflict with 
observations. A quasar could be in the very high state at 
the present time and show no signs of limit cycle behavior, 
as the time scale for the cycle should be 10 s times longer 
than the stellar-mass black hole case (several thousand 
years). Only every few thousand years would a quasar en- 
ter a very active phase that lasted a few decades. During 



that time a temporarily high accretion rate (m ~ 1) pe- 
riod may occur, with jet powers also given by equations 
(H) and (0), while during the quiescent phase the jet pow- 
ers would be given by equations (^) and (||). The final jet 
power on the kiloparsec scale (integrated over thousands of 
years and several limit cycles) would be a sum of the two 
pertinent equations, weighted by the duty cycle of each 
phase — that is, approximately equation (^|) or (^), with 
rh replaced by the accretion rate averaged over the long 
term. 

Accreting supermassive black holes, then, should have 
only two basic states for each of the Schwarzschild and 
Kerr cases: low/hard (ADAF) or very high (geometrically 
thin, radiation-pressure-dominated and unstable flows). 
We consider the very supermassive (10 7_9 MQ/quasar) 
case first. The Schwarzschild low/hard case is described 
only by equation (^J), and corresponds to weak radio 
cores with little or no line e mission ("Class B" accord- 
ing to Jackson & Wall 199E). As discussed above, the 



Schwarzschild very high disk case is a linear combination 
of equations (^) and (p|), yielding a long-term jet power 
of Lf£ w - 5 x 10 43 ergs" 1 f 2 for an ADAF duty cy- 
cle of 1% and a 10 9 A/q hole. In the latter state the hole 
usually would produce significantly more optical-UV ra- 
diation than in the low/hard state, appearing as a radio 
quiet "Class A" object (QSO). The Kerr low/hard case is 
described by equati on ([7|), and correspond s to a Class B 
radio galaxy (either [Fanaroff fc Riley 1974 [FR] type I or 
type II without strong nuclear emission lines). The Kerr 
very high state is again a linear combination of equations 
(§) and (@), yielding Lf e f r ~ 10 45 ergs -1 g 2 for a 1O 9 M 
hole and a 1% duty cycle — a radio loud Class A object 
(quasar or FR II BLRG) with a jet perhaps hundreds of 
times more powerful than in the above radio quiet case. 

In the moderately supermassive (1O 5_7 M0/Seyfert) 
case, significant optical emission should be present in both 
the low/hard and very high states, as that emission would 
be produced in the middle regions of the disk, which are 
not affected by state transitions. In the very high state, 
during the thin-disk portion of the limit cycle, the peak 
emission would be in the extreme ultraviolet (~ 0.1 keV). 
Such objects should appear similar to the highly- variable , 
ultrasoft narrow- line Seyfert 1 galaxies ( Brandt 2000| ), 
whereas in the low/hard emission state they should ap- 
pear more like the classical hard X-ray selected Seyferts. 
Their radio jet properties should be similar to those of 
the quasars above, though 100 times weaker, with the ul- 
trasoft objects producing a jet only when they enter the 
thick-disk portion of their limit cycle. 
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